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Abstract 

We consider the 0{n) loop model on tetravalent maps and show how to 
rephrase it into a model of bipartite maps without loops. This follows from a 
combinatorial decomposition that consists in cutting the 0{n) model configu- 
rations along their loops so that each elementary piece is a map that may have 
arbitrary even face degrees. In the induced statistics, these maps are drawn 
according to a Boltzmann distribution whose parameters (the face weights) are 
determined by a fixed point condition. In particular, we show that the dense 
and dilute critical points of the 0{n) model correspond to bipartite maps with 
large faces (i.e. whose degree distribution has a fat tail). The re-expression 
of the fixed point condition in terms of linear integral equations allows us 
to explore the phase diagram of the model. In particular, we determine this 
phase diagram exactly for the simplest version of the model where the loops 
are "rigid" . Several generalizations of the model are discussed. 



1. Introduction and main results 

1.1. General introduction 

Planar maps, which are proper embeddings of graphs in the two-dimensional sphere, 
are fundamental mathematical objects whose combinatorics raises many beautiful enu- 
meration problems, first addressed by Tutte in the 60's f^. Maps are also widely used 
in physics as discrete models for fiuctuating surfaces or interfaces in various contexts, 
ranging from soft matter physics (e.g. biological membranes) to high energy physics 
(e.g. string theory). Most problems involve "random maps", i.e. consist in the study of 
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some particular statistical ensemble of maps, distributed according to some prescribed 
law. Interesting scaling limits may be reached when one considers ensembles of large 
random maps, giving rise to nice universal probabilistic objects. Several map enumera- 
tion techniques were developed over the years, which we may classify in three categories: 
Tutte's original recursive decomposition method P], the technique of matrix integrals ^ 
and, more recently, the bijective method |^] where maps are coded by tree-like objects. 

So far, the most advanced results were obtained for ensembles of maps with a 
simple control on the degree of, say the faces of the maps. Examples are ensembles 
of triangulations (maps with faces of degree 3 only) or of quadr angulations (degree 
4) where the total number of faces in the map is either fixed or governed by some 
Boltzmann distribution. The statistics of these maps is now well understood and many 
exact enumeration results were obtained within the framework of each of the three 
above enumeration techniques. Of particular interest is the scaling limit of large maps 
with prescribed bounded degrees, which realizes the so-called universality class of "pure 
gravity" , and gives rise to the "Brownian map" , a probabilistic object with remarkable 
metric properties [Q. As first recognized in physics, other universality classes of maps 
may be reached upon equipping the maps with statistical models, such as models of 
spins or particles, which present a large variety of critical phenomena that modify 
the statistical properties of the underlying maps. This gives rise to a large variety 
of universality classes for maps, whose understanding is the domain of the so-called 
"two-dimensional quantum gravity" |^ . It is worth mentioning that most results in this 
domain rely on matrix integral formulations of the models, so that the matrix integral 
technique appears so far as the most powerful approach to explore new universality 
classes. 

A particular important class of models of maps equipped with statistical models 
are the so-called 0(n) loop models which consist in having maps endowed with con- 
figurations of closed self- and mutually-avoiding loops drawn on their edges, each loop 
receiving the weight n. A particular class of 0{n) loop models where loops visit only 
vertices of degree 3 (the degree of the other vertices being bounded) was analyzed in 
details by use of matrix integral techniques in [0-0. There it was found that these 
models present several phases with non-trivial universal scaling limits. We will briefiy 
recall the results of this analysis in Section 1.2 below. 



It was recognized recently in that another simple way to escape from the 
universality class of pure gravity consists in considering "maps with large faces", i.e. 
ensemble of maps where the degree of the faces is unbounded and properly controlled 
so that faces with arbitrarily large degree persist in the scaling limit. Using now the 
bijective method, it was shown that such ensembles may give rise to new probabilistic 
objects corresponding to maps coded by stable trees. It was then proposed that these 
probabilistic objects may also describe the scaling limit of 0{n) loop models in some of 
their phases. 

The purpose of this paper is to show that this is indeed the case. More precisely, we 
show that a number of 0{n) loop models may be bijectively transformed into models of 
maps without loops and with a simple control on the degree of their faces. This control 
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involves some degree-dependent face weights whose value is fixed by some appropriate 
consistency relation in the form of a fixed point condition depending on n. This implies 
that the possible scaling limits of the 0(n) loop models necessarily match the scaling 
limits of maps controlled simply by their face degrees. These include the pure gravity 
universality class as well as the new universality classes of Ref. for maps with large 
faces, which are observed in the dense and dilute phases of the 0(n) loop models. 

1.2. Previously solved 0{n) loop models 

As a reminder, we briefly describe a few of the known results about 0{n) models on 
random maps, as previously studied in [0-0 by matrix integral techniques. We will not 
give very precise definitions here as these models are not those that we will study in the 
remainder of the paper. The 0{n) loop model studied in [0-|l2| consists in having self- 
and mutually-avoiding loops drawn on random maps so as to visit only vertices of degree 
3, the vertices not visited by loops having arbitrary (but bounded) degrees at least 3. 
A vertex visited by a loop receives the weight h while a vertex not visited by a loop 
receives the weight gk if it is /c-valent (/c > 3). Each loop receives in addition the weight 
n. We define -Fp°°P as the corresponding generating function for loop configurations on 
rooted maps, where the root vertex is not visited and has degree p (by convention, the 
root vertex is unweighted and we set Fq°°^ = !)• This corresponds to having a dual map 
with a boundary of length p. We may gather the family (Fp°°P)p>o into the resolvent 

/'loop 

^(0 = E7?Tr' (1-1) 

well-defined for ^ large enough. In the range of weights where the model is well-defined, 
this formal series is in fact an analytic function on C \ [7_ , 7_|_] , and has a discontinuity 
on some segment [7_ , 7_|_] . Using this information in the recursive relation for Fp°°P 
obtained (in the spirit of the Tutte's recursive method) by removing the root edge, it 
can be shown that W{^) is solution of a scalar non-local Riemann-Hilbert problem: 

e [7-, 7+], Wi^ + iO) + - iO) + nWCh-' - = ^ - J] ^fc e'"' • (1-2) 

fc>3 

This fixes uniquely W{^) as well as 7_ and 7-(-, from the requirement that W{^) ~ 1/^ 
when ^ —7- oo, and that W is holomorphic in C \ [7_, 7-(-]. A non-trivial critical point is 
reached when 7_|_ — ?■ (2/i)~^, which defines a non-trivial critical surface h = h{{gk)k>i)- 
The equation (|1.2|) was first solved on this critical surface by Rostov who showed 
that W may develop a singularity of the form: 



Vr(0|sing. ~ const, for , (1.3) 
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with a = 2 ±b, ixh = arccos(n/2). This singularity captures the large p asymptotics of 
Fj^°°P as 

F^°°P~ const, ^^fc. (1.4) 

The smallest value a = 2 — 6 is observed inside the non-trivial critical surface and 
describes the so-called dense phase of the 0{n) model, while the largest value a = 2-\-h 
is observed only at some boundary of this critical surface and describes the so-called 
dilute phase. Other multicritical points, with a = 2 ± 6 + 2m (m = 1, 2, . . .), can also be 
observed by tuning more and more coefficients gk 0, but this requires having some of 
these weights negative, which prevents interpreting these multicritical points as proper 
probabilistic ensembles. Eq. ( p..2|) was then solved in all generality (outside of the critical 
surface) in [jll| in terms of elliptic functions and of the Jacobi theta function -{?(•, g). 



the non-trivial critical surface corresponding to the limit where the nome q goes to 1 
(or depending on the convention). In this limit, the theta function degenerates into 
a trigonometric function, and one can recover the critical exponents stated above. The 



study of the 0{n) model was extended recently to maps of any topology in Ref. [jl4 



1.3. Overview of the paper 

Let us start by giving a brief description of our results. Most of the paper deals 
with 0{n) loop models defined on planar tetravalent maps (maps with vertices of degree 
4 only). By definition, the loops are self- and mutually-avoiding and each loop receives 
the weight n. Different models are obtained by assigning different weights to the vertices 
depending on whether they are visited by a loop or not. A special attention will be paid 
to the simplest version of the model, the rigid loop model, where loop turns are forbidden 
at the tetravalent vertices. 

Our analysis is simplified by first reformulating our loop model on the dual maps, 
which are planar quadr angulations, and by then extending them to quadrangulations 
with a boundary of arbitrary (necessarily even) length. Indeed, a simple characterization 
of the universality class of the loop model at hand is via the large p asymptotics of its 
generating (or partition) function -F^^^p in the presence of a boundary of length 2p. It 
is expected to take the form 

F^°°P ~ const.— (1.5) 

with a non-universal exponential growth factor A and a power law decay factor 1/p" 
involving a universal exponent a in the range 3/2 < a < 5/2. Assuming this asymptotic 
behavior, we show that, for a given n, at most four values of the exponent a may be 
observed, namely a = 3/2, corresponding to a suhcritical vaodel, a = 5/2, corresponding 
to a generic critical model, and 

a = 2 ± 6, 7r6 = arccos(n/2) (1-6) 

corresponding to a non-generic critical model. This last behavior requires that n lies 
in the range < n < 2 (note that, for n = 0, ( p..6|) yields the pure gravity exponents 
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a = 3/2 and a = 5/2). In the case of the rigid loop model, we solve the model exactly, 
thus establishing the validity of ( p. .51) and providing a precise description of the phase 
diagram which specifies the domain of physical parameters where each of the above 
values of a is observed. 



The expressions ( p..5|) and ( |1.6|) are identical to those obtained for the 0{n) loop 



models discussed in Section 1.2. Here, they are obtained by different means which rely 
on an equivalence of our models with models of bipartite maps. More precisely, quadran- 
gulations (with a boundary) endowed with loop configurations may be coded bijectively 
by bipartite maps without loops but with faces of arbitrary degree, each containing an 
0{n) loop model configuration of its own. This bijective decomposition allows, via a 
simple substitution procedure, to express Fp°°P in terms of the well-understood gener- 
ating function for bipartite maps with a boundary. As a consequence, we may identify 
the possible asymptotics of as those of ordinary maps with possibly large faces. 



leading eventually to ( p..5|) and (|1.6|) . 

The paper is organized as follows. Section 2 is devoted to a combinatorial study of 
the 0{n) loop model on tetravalent maps. After defining the model in Section 2.1, we 
present in Section 2.2 a bijective decomposition which allows to code the configurations 
of our model in terms of their gasket, which is a bipartite map with no loops but with 
holes, which are faces of arbitrary even degree, together with a content for each hole. 
For consistency, the effective weights for the holes are shown to obey some crucial fixed 
point condition which determines them uniquely. This condition involves in particular 
some ring generating function, which accounts for the configuration in the immediate 
vicinity of a loop, and which we make explicit in Section 2.3. Section 3 does not treat 
the 0{n) model itself, but discusses the possible critical behaviors of bipartite maps, 
which we classify into subcritical, generic critical and non-generic critical. The first two 
cases are discussed in Section 3.2 while a particular attention is paid to the non-generic 
critical behavior in Section 3.3. Section 3.4 discusses the so-called resolvent of bipartite 
maps. Section 4 is devoted to the consequences of the fixed point condition on the 
asymptotics of Fp""^, leading to ( P^ ) and ( |1.6| ) above. We first concentrate on the rigid 
loop model in Section 4.1 before addressing the general case in Section 4.2. Section 5 
shows how to transform the fixed point condition into a linear integral equation, both 
in the rigid case (Section 5.1) and in the non-rigid case (Section 5.2). We then analyze 
in Section 5.3 the solution of this equation as it captures all the phase diagram of the 
0{n) loop model. Section 6 presents a detailed analysis of the rigid loop model. We 
first show in Section 6.1 how to transform the fixed point condition into an equation for 
the resolvent of the model. This equation is then solved, first along a line of non-generic 
critical points in Section 6.2, then in all generality in Section 6.3. The generic critical 
line is discussed in Section 6.4 and we summarize the resulting phase diagram in Section 
6.5. Section 7 considers extensions of our results to other classes of 0{n) loop models. 
We first show how the relation ( |1.6| ) may be modified by considering loops with non- 
symmetric weights (Section 7.1) or with additional constraints on their length (Section 
7.2). We end our study by a discussion of the wide class of 0{n) loop models defined 
on arbitrary even-valent maps (with bounded degrees) and show that our results nicely 
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extend to this case. We gather our conclusions in Section 8. 



2. The 0{n) loop model on tetravalent maps: combinatorics 

2.1. Definition of the model 

Our model of interest is the 0{n) loop model on planar tetravalent maps defined 
as follows. By loop we mean an undirected simple closed path on the map (visiting 
edges and vertices), also sometimes called an undirected cycle. A loop configuration 
is a set of disjoint loops. By construction, loops are both self- and mutually-avoiding. 
Alternatively, a loop configuration may be viewed as a set of covered edges such that 
each vertex is incident to either or 2 covered edges. 





(a) 



(b) 




Fig. 1: A configuration of the 0{n) loop model on tetravalent maps, when 
viewed on the dual quadrangulation, is built out of three types of squares: 
(a) empty squares, weighted by g, (b) squares visited by a loop going 
straight, weighted by hi, and (c) squares visited by a loop making a turn, 
weighted by /i2- 



We prefer to work with the dual map which is a planar quadrangulation. Then, the 
loops cross some edges of the quadrangulations. Up to a rotation, the quadrangulation 
has three possible types of faces, as displayed in Fig. 1: (a) empty squares not visited 
by a loop, (b) squares visited by a loop going straight (i.e. crossing opposite sides of 
the square) and (c) squares visited by a loop making a turn (i.e. crossing consecutive 
sides of the square). The length of a loop, i.e. the total number of faces that it visits, 
may be even or odd. However, the number of faces of type (b) visited by a given loop 
is necessarily even, resulting from the fact that a planar quadrangulation is bipartite 
(other interesting consequences of the bipartite nature of our maps will be seen below). 
A loop is said rigid if it only visits faces of type (b). 

Generally speaking, the 0{n) loop model consists in attaching a non-local weight 
n to each loop, in addition to some local weights. Here, we consider an "annealed" 
model where the map and the loops are drawn at random altogether: a configuration of 
the model is the data of a planar tetravalent map endowed with a loop configuration. 
Local weights are attached to the vertices of the tetravalent map, or equivalently to the 
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faces of the dual quadrangulation, and naturally these weights depend on the face type, 
see again Fig. 1: g per square of type (a), hi per square of type (b), /i2 per square 
of type (c). The global weight of a configuration is the product of all loop and face 
weights, and the generating ( "partition" ) function of the model is the sum of the global 
weights over all configurations. Here n, g, hi, /i2 are taken as non-negative real numbers. 
The model is said well-defined when the generating function is finite (this is the case 
for instance when g + max(n, l)(2/ii + 4/12) < 1/12). It is then possible to normalize 
the configuration weights and define a probability distribution. The rigid loop model 
corresponds to taking /i2 = so that loops are obliged to go straight within each visited 
square. 




Fig. 2: A loop configuration on a quadrangulation with a boundary of 
length 2p = 8, with 3 loops, 9 squares of type (a), 8 squares of type (b) 
and 11 squares of type (c). Its weight is thus n'^g'^h\h^. The associated 
gasket consists of 5 regular squares and 2 holes with degrees 6 and 10. 

As usual it is easier to work with rooted maps. Recall that a rooted map has a 
distinguished oriented edge (the root edge). From now on, the face on the right of the 
root edge will be called the external face and we define a quadrangulation with a boundary 
of length 2p (p > 1) as a rooted bipartite planar map where the external face has degree 
2p and every other face has degree 4. We shall now consider loop configurations on the 
dual of a quadrangulation with a boundary (see Fig. 2) and, for simplicity, we assume 
that no loop visits the external face, to which we therefore decide to attach a local 
weight 1. Keeping the same weights n,g,hi,h2 as above, we define Fp°°P(n; (7, /ii, /i2) 
as the generating function for the 0{n) loop model on duals of quadrangulations with 
a boundary of length 2p. As particular cases, for p = 1 we obtain (upon squeezing the 
bivalent external face) rooted quadrangulations where the root edge is not crossed by a 
loop, while for p = 2 we obtain (directly) rooted quadrangulations where the external 
face is not visited. 
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2.2. The gasket decomposition 

The fundamental observation of this paper is that F^°''P{n]g, hi,h2) is related the 
generating function Fp of bipartite planar maps without loops with a boundary of length 
2p, defined as follows. Let us recall that a planar map is bipartite if and only if all its 
faces have even degree. We say that it has a boundary of length 2p if it is rooted and 
the external face, defined again as the face on the right of the root edge, has degree 2p. 
To each internal face of degree 2/c, (/c > 1), we attach a weight g^, and the weight of a 
rooted map is the product of all its internal face weights. We then define Fp as the sum 
of the weights of all bipartite maps with a boundary of length 2p. Our main statement 
is that 

Fl--^{n-gMM) = Fp{g^,g2....) (2.1) 
for the particular sequence {gk)k>i of face weights satisfying the fixed point condition 

gk = gSk,2 + n^ Ak,k'{hi,h2)Fk'{gi,g2,---), k>l, (2.2) 

fc'>0 

with Ak^k'i^i, /i2) a polynomial in hi and /i2 whose expression will be given in the next 
Section. For the rigid case, it reads simply Ak,k'{hi, 0) = h'l''Sk,k', leading to the simpler 
fixed point condition 

9k = gh,2 + nhl^Fk{gi,g2-, • ■ •) • (2-3) 




Fig. 3: The outer and inner contours of a loop. 

Let us now justify our fundamental observation, by considering a quadr angulation 
with a boundary of length 2p endowed with a loop configuration. For conciseness, 
we will often, in this Section, call simply "quadrangulation" a quadrangulation with a 
boundary. Thanks to the external face, the notions of exterior and interior of a loop 
are well-defined. To each loop, we may associate its outer (resp. inner) contour formed 
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by the edges of the quadr angulation which (i) belong to the exterior (resp. interior) 
of the loop and (ii) are incident to squares visited by the loop (see Fig. 3). Note that 
these contours are closed paths living on the quadrangulation itself and therefore each 
of them has an even length. Note also that the inner contour may be empty, i.e. may 
have length if the loop encircles a single vertex. 




Fig. 4: The gasket is obtained by wiping out the content of the outer 
contours of the outermost loops, creating holes. 



Using a terminology similar to |1T3|] , we call gasket the map formed by edges that 
are exterior to all the loops, see Fig. 4 for an illustration. It is a rooted bipartite planar 
map containing in particular all the edges incident to the external face, which therefore 
remains of degree 2p. Its other faces are of two types: regular faces which are the squares 
of the original quadrangulation that lay outside of every loop, and holes delimited by 
the outer contours of the former outermost loops. 

Clearly, the transformation that goes from the original quadrangulation with loops 
to the gasket is not reversible. To make it reversible, we must keep track of the former 
content of the holes. In this respect, a first remark is that the interior of a loop may itself 
be viewed as a quadrangulation with a boundary endowed with a loop configuration, 
where the boundary is nothing but the inner contour of the loop at hand (when the 
inner contour has length 0, we have the so-called "vertex-map" reduced to one vertex 
and one face, with no edge) . In order not to lose any information, we must also consider 
the squares covered by the loop itself, which form a ring lying in-between the inner 
and outer contours. As shown in Fig. 3, this ring is a cyclic sequence of squares glued 
together. To summarize, the content of a hole of degree 2k {k > 1) is described by a 
pair formed by a ring with outer length 2k and inner length 2k' (for some k' > 0), and 
an internal quadrangulation with a boundary of length 2k' (equal to the vertex-map if 
k' = 0) endowed with a loop configuration. 

Two important remarks are in order. First, to ensure bijectivity, both the ring 
and the internal quadrangulation must be rooted, respectively on its outer contour 
and on its boundary. The positions of their root edges is inherited from the rooting 
of the original quadrangulation via a somewhat irrelevant yet well-defined procedure: 
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Fig. 5: Illustration of the rooting procedure. We draw the leftmost short- 
est path (dashed line) in the gasket starting with the map root edge and 
ending with an edge of the outer contour of the loop at hand, which we 
pick as root for the ring. We then select an edge of the inner contour by the 
rules displayed at the bottom left, which we pick as root for the internal 
quadrangulation. 

for instance we may consider the leftmost shortest path that stays within the gasket, 
starts with the root edge and ends with an edge of the outer contour, that edge being 
selected as the ring root; then, using only the data of the ring, we may easily select 
an edge of the inner contour which serves as root for the internal quadrangulation 
(see Fig. 5). Second, an essential assumption is that the original map, the gasket and 
the internal quadrangulations are all possibly separable (i.e. may contain separating 
vertices whose removal disconnects them). In particular, a separating vertex incident 
to a hole in the gasket appears whenever a multiple point was present along the outer 
contour of the associated loop on the original quadrangulation. Similarly, a multiple 
point of the inner contour of a loop results into a separating vertex on the boundary 
of the associated internal quadrangulation. In contrast, all the vertices of a ring are 
considered as distinct (since the information about contacts is recorded in the gasket 
and internal quadrangulation). At this stage, it should be clear that the decomposition 
is reversible: given a ring and a quadrangulation with a boundary with compatible 
lengths, there is a well-defined procedure to join them and fill a hole of the gasket. 

Globally, we have a one-to-one correspondence between, on the one hand, quadran- 
gulations with a boundary endowed with a loop configuration and, on the other hand, 
gaskets endowed with hole contents. Translating this correspondence in the language 
of generating functions, we deduce that the wanted quantity Fp°°P(?T,; (7, hi, is equal 
to the generating function Fp for gaskets, where each regular face still receives a weight 
g and each hole of degree 2k receives a weight nY^Ak^k'{hi,h2)F]^,°^{n;g,hi,h2) (see 
Fig. 6). In this hole weight, n accounts for the associated loop, F^°^{n] g, hi, /i2) for the 
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Fig. 6: A face of the gasket is either a regular square (weighted by g) 
or a hole of even degree, say 2k. The content of this hole is made of a 
ring with outer length 2k and inner length 2k' for some k' > (weighted 
by nAk^k'i^i, h2)), and of an internal quadr angulation with a boundary of 

length 2k' (weighted by F^^°^{n] g, hi, /i2) = Fk>). 



internal quadrangulation, while Ak^k'{hi,h2) is the generating function for rings with 
sides of lengths 2k and 2k' . This establishes Eqs. ( pTlD and (pT2D. 

2.3. The ring generating function 

Let us now discuss the precise form of the ring generating function ^/(/ii, /i2)- 
We have the explicit expression 

-M)^ E (2.4, 

which follows from a simple counting argument: the term with index j in the sum 
corresponds to rings with 2j squares of type (b), each weighted by hi. Such rings have 
k — j (resp. k' — j) squares of type (c) with their two unvisited sides along the outer 
(resp. inner) contour of the ring, all having a weight /i2- There are k-j^k'-j) possible 
sequences of such three types of squares. However, this number has to be corrected 
because a ring has a distinguished edge (among 2k) on its outer contour rather than a 
distinguished square (among k + k') to start the sequence. This explains the corrective 
factor 2k/{k + k'). 

We have the alternative expression, which will prove useful in the following: 
E A,,,,(h„h,)z-''MM(z)f HX-(z)f , A±(.)= . ± Vfej-4/.i+4ft. .y 

A;'>0 \ \Z 2) J 

(2.5) 

(note that the sum converges for \z\ > h2). This expression naturally follows from a 
transfer matrix argument. As seen in Fig. 7, we may distinguish two "states" for the 
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Fig. 7: The ring transfer matrix M{z). An edge of the outer contour 
(thick edges) is in state 1 if it follows clockwise a vertex adjacent to a 
vertex of the inner contour, and in state 2 otherwise. The dashed lines 
indicate how we distribute the weights in the transfer matrix. Each edge 
in state 2 is followed by an edge in state 1 and receives the weight /i2 of 
its incident square of type (c). Each edge in state 1 first receives a weight 
1/(1 — h2Z~^) accounting for a sequence of squares of type (c) as shown 
here in grey, then a further weight hiz~^^'^ if and only if it is followed by 
another edge in state 1. 

edges of the outer contour depending on whether they follow clockwise a vertex adjacent 
to a vertex of the inner contour (state 1) or not (state 2). We may then write 



Afc,fc,(/^l,/^2)^-^'=trM(;^)2^ M{z) = 



k'>0 



l-h2Z-^ l-h2Z-^ 
\ h2 / 



(2.6) 



where M(z) is the transfer matrix describing the transition between two successive edge 
states. This establishes ( |2.5| ) where X±{z) are simply the eigenvalues of M'^{z), solutions 
of: 

{h2{X + z) - Xzf - hlXz = . (2.7) 

Interestingly, A and z play a symmetric role in this equation. It follows that A(A(z)) = z 
provided we pick the correct determination of A. More precisely, it may be checked that 
z I—)- A_|_(z) is an involution on (/125 00) with a unique fixed point: 

z* = hi + 2h2 = X+{z*) . (2.8) 

Finally, we observe that for rigid loops, there is a unique ring with outer length 
2/c, and it has the same inner length. This yields immediately ^/(/ii, 0) = hi^dk,k' as 
announced in the previous Section. This is compatible with (|2.4|) for h2 = 0, and with 
( ^75|) upon noting that X+{z) = h\/ z and X-{z) = 0. 
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3. Bipartite maps with arbitrarily large faces: possible critical behaviors 



In this Section, we step away from the 0{n) loop model and discuss a few properties 
of the generating functions of bipartite maps. 

3.1. Reminders 

As already mentioned, a planar map is bipartite if and only if all its faces have even 
degree. To each face of degree 2/c, (/c > 1), we attach a weight gk, here considered as 
an arbitrary non-negative real number. Furthermore, we will sometimes find convenient 
to also attach a non-negative real weight u per vertex of the map, irrespectively of its 
degree. The unnormalized weight of a map is then the product of all its face and vertex 
weights. Note that, from Euler's relation, we could set u = 1 without loss of generality 
upon redefining gk gku'^~^ ■ 

It is well-known that the simplest map enumeration formulas are for pointed rooted 
maps, i.e. maps with a distinguished vertex (pointed) and a distinguished oriented edge 
(rooted). The generating function (i.e. the sum over all maps of the unnormalized 
weights above) for pointed rooted bipartite planar maps is 2u{R{u) — u) where R{u) is 
the smallest non-negative (possibly infinite) solution of the equation Jir 



fc>i 

Here, we choose to display explicitly the dependence in u but to hide that in the gus. 
It is easily seen that R{u) is an increasing function of u with -R(O) = 0. 

Other more involved generating functions can be expressed in terms of R{u). For 
instance, the generating function for maps with a boundary of length 2/c, as defined 
in the previous Section (with u = 1), reads 



2k 
k 



1 

k 



R{uYdu . (3.2) 







Consistently, we set Fq = 1. Furthermore, since rooted maps are obviously in bijection 
with maps with a boundary of length 2 and not reduced to a single edge, we may 
interpret Fi — 1 as the generating function for rooted maps. Upon integrating (|3.1| ) and 
noting that {^^^ = 2(^'^^^), we immediately deduce 

F^-l = Y,9kFk (3.3) 

fc>i 

which amounts to decomposing the generating function of rooted maps (l.h.s) according 
to the root degree of the maps (r.h.s). 
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3.2. Subcritical and generic critical ensembles 

We say that the sequence {gk)k>i of face weights is admissible if -R(l) is finite. In 
this case, taking u = 1 and dividing the unnormahzed weight of a pointed rooted map 
by 2(i?(l) — 1), we obtain its probabihty in the Boltzmann ensemble of pointed rooted 
maps associated with the sequence {gk)k>i- For example, random quadrangulations are 
obtained by choosing the sequence gk — gSk,2 (fi* > 0), which is admissible for g < 1/12. 

A number of properties of the Boltzmann ensemble are encoded in the expansion 
of R{u) around u = 1. This expansion is obtained upon inverting the expansion around 
R = -R(l) of the inverse function u{R) which, from ( |3.1| ), reads explicitly: 

u{R) = R-^{R), ^{R) = Y,^gk(^^~'^^R' . (3.4) 

Here the function f>{R) is defined through its series expansion at i? = 0, whose coeffi- 
cients are all non-negative. Its radius of convergence is Rc = 1/(4 limsupj[.^^((7fc)^/'^) 
since (^^^) ~ 4'^/v^7tA; at large k. If the sequence is admissible, Rc is necessarily non- 
zero (possibly infinite), with -R(l) < -Rc, and it is easily seen that u'{R{l)) > 0, i.e. 
(p'{R{l)) < 1. The sequence is said critical if (p'{R{l)) = 1 (hence u'{R{l)) = and 
R'{u) — i- oo as w — )■ 1) and subcritical otherwise (hence -R'(l) is finite). This change 
of behavior is visible in the large k asymptotics of F^. Indeed, the integral (|3.2| ) is 
asymptotically dominated by the vicinity of m = 1 and its behavior may be computed 
via Laplace's method |]T^. For a subcritical sequence, we have at large k: 

p_ R{i){i-^'{R{i)){m{i)f 



This form of the asymptotics is observed for instance when all gkS are zero, in which 
case 9?(-R) = 0, R{u) = u, and Fk = [^^) /{k + 1) is the k-th Catalan number. 

For a critical sequence, the asymptotic behavior depends on whether the (negative) 
quantity u"{R{l)) = —<f>"{R{l)) is finite or not. For a generic critical sequence, this 
quantity is finite and we have that 

R{1)V{R{1)) {4R{l)r 
^ k^ ■ ^'-'^ 



This form of the asymptotics is observed for instance for critical random quadrangula- 
tions obtained by choosing the sequence gk = (1/12)5^^2, in which case (p{R) = R^/4:, 
R{1) = 2, Lp"{R{l)) = 1/2, so that ( p.6| ) is consistent with the exact formula 

I? _ ofc+l (2fc)! 
— ^ fc!(fc+2)!- 

In order to escape from this generic critical behavior, we must have (/?"(i?(l)) in- 
finite, which implies that -R(l) must coincide with the radius of convergence Rc of (f. 
Such non-generic critical behaviors are analyzed in the next section. Before proceeding, 
let us mention the following simple lemma: 
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• Monotonicity property: if {gk)k>i and {gk)k>i are two sequences such that gk < (jk 
for aU k, the inequahty being strict for at least one k, and if {gk)k>i is admissible, then 
{gk)k>i is subcritical. 

Indeed, denoting by R{-) and (^(■) the functions associated with the sequence 
{gk)k>i, we have -R(l) < -R(l) since 2R{1) — 1, being a map generating function, is 
a strictly increasing function of gk for each k. Then, (f>'{R{l)) < (f>'{R{l)) < 1 hence, 
by definition, {gk)k>i is subcritical. 

3.3. Non-generic critical ensembles 

In order to get a non-generic critical behavior, we must ensure simultaneously the 
two conditions f'{R{l)) = 1 (criticality) and ip"{R{l)) = oo (non-genericity) . This 
turns out to highly constrain the possible form of the sequence (5'fc)fc>i. As shown by 
Le Gall and Miermont [^, a natural candidate for such a sequence is: 

^ \ k-l 



where {gl)k>i is an arbitrary sequence of reference, such that 

3 5 

gl ~ for /c cx) with - < a < - . (3.8) 

Note that Rc is indeed the radius of convergence of if{R) since we have if{R) = 
cRfoiR/{4:Rc)), where 



fc>i 



is an analytic function with radius of convergence 1/4. Since a lies in the range ]3/2, 5/2[, 
we have (p'{Rc) < oc and (p"{Rc) — oo. The non-genericity condition therefore reduces 
to demanding that 

R{1) = Rc , (3.10) 

or equivalently 1 = Rc — (p{Rc)- Together with the criticality condition ip'[Rc) = 1, this 
fixes the values of c and Rc as 

fl,.l + lMVl). (3.11) 



4/o(l/4) + /i(l/4)' ' /i(l/4) 

At these values, by standard transfer theorems [|18[, u{R) has a singular expansion 
around R = Rc given by 

„(.H:-H-Ii^(.-|)^-'\0(l-|)\ (3.1. 
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Note that r(l/2 — a) > since 3/2 < a < 5/2. Inverting this expansion, we see that 
the generating function R{u) behaves, when we fix the face weights gk as above and let 
the vertex weight u tend to 1 from below (i.e. we use the parameter u to control the 
approach to the critical point), as 

«(„) = «.(i-.(i-»)V<-..)+o(i-.), «-( ,,^^r;f;,_„) )^"^°"' . (3.13) 

Again, by transfer theorems, this implies that the coefficient of in R{u), hence the 
probability of having vertices in the Boltzmann ensemble of pointed rooted maps, 
decays as A^-(2a+i)/(2a-i) as A ^ oo. 

Returning to the asymptotics of Ffc, we find by Laplace's method that, for large k: 

. ^. , 2ci?e(a-l/2)r(l/2-a)r(a-l/2) i^g 
R{u) du (3.14) 

and, using {^^) ~ i^/V^, (a - l/2)r(l/2 - a)T{a - 1/2) = 7r/sin7r(a - 3/2) and 
Rc^R{l), 

2cRil) {4R{1))^ 
smn{a-3/2) k'^ ' ^^'^^^ 

It is instructive to compare this expression with the "input" , namely the form (|3.7] )-( p^) 
for non-generic critical face weights, which implies the asymptotics 

g,^4ci?(l) ^^^^y ■ (3.16) 

We note a striking similarity: the prefactor is simply divided by 2sin7r(a — 3/2) while 
the exponential factor is inverted. This similarity is not merely anecdotal but will have 
important consequences for the 0{n) loop model (see Section 4). 

As a concluding remark, note that the condition i?(l) = Rc for non-generic critical- 
ity implies that the face degree distribution has a fat tail. As seen from the derivation 
|T^ of (|37l]), the probability that the external face in a rooted pointed map has degree 
k is 

R{\) - 1 

which decays exponentially as {R{l)/Rc)^ when R(l) < R^ but only as a power law 
^-a-i/2 ^ non-generic critical point. Similarly, from (|3.3| ), the probability that the 
external face in a rooted (but unpointed) map has degree k is 

_9kFk_ ^3_^g^ 

which also decays exponentially as {R{l)/Rc)^ when < Rc and as a power law, 

now /t , at a non-generic critical point. Note finally that the condition -R(l) = Rc is 
necessary but not sufficient for non-generic criticality as one may construct examples 
of subcritical and generic critical sequences such that -R(l) = Rc- Such situations will 
however not appear in the context of the 0{n) model with non-negative weights. 
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3.4- The resolvent 

For the purposes of Section 6, it is useful to gather facts about the so-called resolvent 
(this terminology being borrowed from the matrix integral formalism), defined as: 

w(o=j:^- (3-19) 

fc>0 ^ 



Using ( |3.2| ) and performing the change of variable u ^ R, we obtain 



WfO^r^H^j. (3.20) 

Wo V 1 - 4i?/^2 

Prom this expression, it is seen that W{^) is analytic in the complex plane minus the 
segment [—7,7], with 7 = 2a/-R(1). Along this segment, it has a discontinuity encoded 
into the so-called spectral density 

2i7r TT J ^2 /4 - ^2 

Since v!{K) > for < -R(l), p(0 is nonnegative on [—7,7]. Note that, by applying 
the Cauchy formula around the cut, we may recover W{£) from p(^) via 

W^(0 = J\y^' ^^[-7,7] (3.22) 

and in particular, since W{^) ~ 1/^ for ^ — )■ 00, we have the normalization J^^ piO^^^ — 
1. The asymptotics (|3.5|), (p.6| ) and ( p.l5|) translate into the respective singularities of 
piO at e = ±7 = ±2y^R{^ 

p(^) ^ ^-^'W)) (4^^i) _ ^2) 1/2 g^bcritical 

271" 

P(0 ~ ^^^7^^ (4i?(l) - ef^^ generic critical (3.23) 

P(0 ~ ^ • ; ^(4^(1) - i^T'^ non-generic critical . 

r(a) sm7r(a — 3/2) 

4. Asymptotic self-consistency 

In this Section, we return to our 0{n) loop model and make a first analysis of the 



fixed point condition (|2.2|) via its asymptotic consequences. The fixed point condition is 
an equation for the unknown sequence {gk)k>i that depends on the physical parameters 
71, g, hi and /i2- By definition, the model is well-defined when a finite solution exists 
(i.e. Fp°°P is finite). Since it involves the generating functions {Fk)k>Q, such a finite 
solution {gk)k>i is necessarily admissible in the sense of Section 2. The purpose of this 
Section is to classify the possible nature (subcritical or critical, generic or not) of the 
solution. 
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4.1. Rigid case 

For simplicity we first concentrate on the rigid case. In this Section, {gk)k>i de- 
notes a sequence satisfying the fixed point condition ( ^.31) . It immediately implies the 
asymptotic condition 

gk ~ n/if Ffc {k oo). (4.1) 

In Section 3, we have classified the possible asymptotic behaviors of Fk depending on the 
nature of the weight sequence {gk)k>i- the equations ( |3.5|) , ( |3.6| ) and ( p.l5|) describe 
respectively the subcritical, generic critical and non-generic critical cases. Here we 
assume that, in the non-generic critical case, the weight sequence is of the form ( p.7|) - 
( p.8|) . It is shown in Section 6 that no other more exotic non-generic behavior exists. 

We now remark that Eq. ( |3.2| ) implies that (Fk)^^'^ — )■ 4i?(l) for large k, regard- 
less of the nature of the weight sequence. By comparing with the general relation 
limsup;j_j.oo(5rfc)^/^ = l/(4i?c), where Rc is the radius of convergence of the function 
ip{R) defined in Eq.(|3.4|), Eq.( [4.1| ) implies the relation 

16hlRcR{l) = 1 . (4.2) 

For later use, we introduce the parameter 

r = 4hiR{l) . (4.3) 

Since -R(l) < R^ we have by (|4.2| ) that r < 1 and the condition r = 1 amounts to 
R{1) = Rc- We now prove the following: 

• Non-genericity criterion: -R(l) = Rc (i-e. r = 1) if and only if the sequence {gk)k>i 
is critical and no n- generic. 

The first implication ("if") was shown in Section 3 without recourse to the fixed 
point condition (see Eq. (|3.10| )). Conversely, assume that the sequence {gk)k>i is either 
subcritical or generic critical. Its asymptotic behavior is directly deduced from that of 
Fk via Eq.([4.1D. In the subcritical case, due to the /c~^/^ factor in Eq.(|3.5|), ip'{Rc) is 
infinite while, by definition of subcriticality, we have <f>'{R{l)) < 1, hence -R(l) 7^ -Re- 
in the generic critical case, due to the /c~^/^ factor in Eq. (|3.6|) , (p"{Rc) is infinite while, 
by definition of genericity, we have (p"{R{l)) < 00, hence -R(l) 7^ Rc- This completes 
the proof (note that when we do not impose the fixed point condition, we can easily 
construct examples of subcritical and generic critical sequences such that -R(l) = Rc)- 
In view of the discussion at the end of Section 3.3, we deduce that the face degree 
distribution of the gasket, i.e. the outer loop length distribution in the 0{n) model, has 
a fat tail iff the model is at a non-generic critical point. 

If we now assume that the sequence {gk)k>i is critical and non-generic, and takes 
the general form ( p77| ) for a suitable sequence {gl)k>i satisfying (|3.8| ), then, in addition 
to the condition r = 1, we see by comparing the prefactors of ( |3.15|) and ( p.l6|) that 

n = 2sin7r(a-3/2) . (4.4) 
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In particular, n must be in the range ]0, 2[ for non-generic criticality to be possible. 

To summarize, we have shown that the fixed point condition (|2.3|) is compatible with 
the model being either subcritical, generic critical, with r < 1 or, more interestingly, 
non-generic critical with r = 1 and a = 2 ±b, nb — arccos(n/2). 

So far, we did not tell which behavior corresponds to a given value of the "physical" 
parameters n, g and hi. Since the generating function Fp'^'^^ {n; g , hi, 0) is an increasing 
function of its parameters, there exists a decreasing function hc{n; g) such that the 
model is well-defined for hi < hc{n;g) and ill-defined for hi > hc{n;g). In particular, 
hc{n;g) > if g < 1/12. By the monotonicity property of Section 3.2, we see that the 
model is subcritical whenever hi < hdji; g). In particular, any critical point necessarily 
lies on the line hi = hc{n; g). An exact expression for hc{n, g) will be given in Section 6 
for < n < 2. There we find that the model is critical all along the line hi = hc{n;g). 
Furthermore, along this line, there exists a g* such that (i) for g* < g < 1/12 the model 
is generic critical, (ii) for < g < g* the model is non-generic critical with a = 2 — b 
{dense model) and (iii) aX g = g* the model is non-generic critical with a = 2 + b {dilute 
model). The corresponding phase diagram is displayed in Fig. 12 below. 

4-2. Non-rigid case 

Let us now consider the general case of arbitrary non-negative values of hi and /i2 
and analyze the consequences of the fixed point condition (|2.2| ) on the asymptotics of gk- 
To treat the subcritical, generic critical and non-generic critical cases simultaneously, 
we write the asymptotics of Fk in the form 

with now 3/2 < a < 5/2 and where x i^ay be read off Eqs.(|3.5|), ( |3.6| ) and ( |3.15| ) 
respectively. Here again, we assume that, in the non-generic critical case, the weight 
sequence is of the form (p.7| ) - (|3.8|) . 

Let us introduce the generating function 

F{z) = Y,Fkz', (4.6) 

fc>0 

related to the resolvent via W{^) = F(l/^^)/^. By transfer theorems, the function F{z) 
is singular when z ^ Zc with a singular part given by 

i^(^)|smg~xr(l-a) (l--) (4.7) 



z. 



where we temporarily exclude the case a = 2. 

To estimate the r.h.s of (^3), we use Eq.(|2.5|) to write 



J2Ak,k'{hi,h2)Fk>= <f^F{z) (A+(z)'= + A_(z)'=) 

k'>0 '' 



(4.8) 
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Fig. 8: Deformation of tiie contour of integration in Eq. (4.8): starting 
from a circle with radius between /i2 and Zc, we deform it by encompassing 
the cut of the function F{z), as shown on the right, so that the dominant 
contribution for /c — > oo is given, after the change of variable z — )■ 5 of 
Eq. (4.9), by integrating the leading order of the integrand over the Hankel 
contour shown below. 

where me may choose as integration contour the circle \z\ = zq for any h2 < zq < z^ 
(note that, since Afc^fc/(/ii, h-i) > Afc^fc/(0, h'l) = , Zc has to be larger than 

/i2 in order for the l.h.s of ( [4.8|) to be finite, hence the model to be well-defined). At 
large /c, this quantity may be evaluated via the method of Hankel contours. Indeed, the 
contour may be deformed as shown in Fig. 8 so that the dominant contribution to the 
integral comes from the vicinity of Zc- Using the scaling 



the variable b is now to be integrated back and forth from cxd to (below and above the 
cut of -F(^)). Using the expansion 




(4.9) 




(4.10) 



with 




(4.11) 



we get at large k 



Ak,k'ihi,h2)Fk' ~ xr(l - a) 



A+(Zc) sin7r(l — a) 



I 

Jo 



oo 



k'>0 



(4.12) 



A+(^c)" 
(;U+)-/c° 
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Here we assumed that X+{zc) > X-{zc), which is vahd as soon as hi > 0. The case hi = 
wiU be discussed later. From the fixed point condition, this behavior is transferred 
into the asymptotics of gk 

^fc~nx^^t^. (4.13) 

By a shght modification of our reasoning, one may easily show that this asymptotics 
also holds in the case a = 2. Using again the property that {gkY^^ — J- l/(4i?c) for 
A; — )■ oo, this yields a first consistency relation 

^ il^^ ■ (4-14) 



This relation generalizes ( [4.2|) which is recovered by noting that \+{z) = h\/z in the 
rigid case. 

Repeating the proof of previous Section, we have here again -R(l) < Rc and -R(l) = 
Rc if and only if the sequence {gk)k>i is non-generic critical. Focusing on such a non- 
generic critical solution, the relation ( |4.14| ) allows to identify l/(4i2(l)) as the fixed 
point z* of the mapping z t-)- A_|_(z), namely 

^ z* = hi+2h2. (4.15) 



4R{1 



A second consistency relation is now obtained by comparing the prefactors in (^4.13|) and 



( p.l6| ) and using the precise value x = 2cR{l)/ sin7r(a — 3/2). We eventually deduce 
the relation 



n = 2(;U*)"sin7r(a-3/2), /U* = - {z^logX+{z) 



-y+iz*) . (4.16) 



Now we recall that z i— )■ X+{z) is an involution in the vicinity of its fixed point z*, which 
implies that X'_^_{z*) = —1 (since X+{z) obviously decreases with z), a result which may 
trivially be checked by a direct calculation. This fixes fj,* = 1 and we recover the same 
relation ( [(.41 ) as in the rigid case. 

To summarize, the fixed point condition ( |2.2| ) is compatible with having the model 
subcritical, generic critical or non-generic critical with a = 2 ± 6, nb = arccos(n/2) 
as long as hi > 0. The case hi = leads to a different relation since in this case 
X+{z) = X-{z) = h2z/{z — /i2) hence the contribution of X-{z) in (f4.8| ) cannot be 
neglected. Then Eq. ( [4.12|) must be replaced by 

E^M'(0,MF.,~2x^^. (4.17) 

At a non-generic point, we get the same relation Zc = l/(4i2(l)) = z* as before (with 
now z* = 2/12) and the same value //* = 1, but comparing the prefactors as was done 
above now yields the new relation 

n = sin7r(a-3/2) , hi = (4.18) 
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without the factor 2 in front. We thus expect that our 0{n) loop model for hi = and 
< n < 1 should have the same non-generic critical behavior as the 0{2n) loop model 
for hi > 0. 



5. Linear integral equation 

In this Section, we show how to rephrase the fixed point conditions (|2.2|) and ( |2.3|) 
as linear integral equations, and then analyze their solutions. We first focus on the rigid 
case before addressing the general case of arbitrary hi and /i2- 

5.1. Derivation of the equation in the rigid case 



We start from the general expression ( |3.4| ) for the function u{R), defining implicitly 
the generating function R{u). Differentiating ( |3^ ) with respect to -R, we readily get 

u'{R) = l-Y,9k(^^~%R'-'- (5.1) 

k>l ^ ^ 

Independently we observe that, upon performing the change of variable u ^ R = R{u) 
in ( 3.2|) , we have 

/ R^u'{R)dR. (5.2) 



So far we have only rewritten some equations of Section 3.2, related to bipartite maps. 
Let us now combine them with the fixed point condition ( p.3|) transcribing the gasket 
decomposition for the 0{n) rigid loop model: assuming that n, g and hi are such that 
the model is well-defined, we plug (|5.2| ) into (|2.3| ), then the result into and obtain 

a linear integral equation for u'{R): 

u'{R)^l-&gR- -^ki\ hf R''-^ I&u'{R)dR (5.3) 

fc>i ^ ^ 



(note that ( ^ ) = (^)/2). By the change of variables 

^=^y y=^y /(^)=«'(^)' P = Q9R{1). r = AhiR{l), (5.4) 
the linear integral equation is rewritten in the more compact form 

fix) + — / Kr{x, y)f{y)dy = 1 - px (0 < x < 1 

^TT Jo 

Krix.y) =T^yij{T^xy), '^(^) = ^Y¥\kl 

k>i ^ ^ 



2 

.fc-1 



(5.5) 
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We recognize a Fredholm integral equation of the second kind \ whose unknown is the 
function / defined on the range [0, 1], and which depends on the non-negative parameters 
n, p, T. Note the consistency relation 



/(-) = ^ (6.6) 

which foUows from 1 = u{R{l)) = /^f u'{R)dR = R{1) f{x)dx. 

Let us now discuss the conditions on the function / and the parameters p and r 
that arise from our derivation. First, the existence of the inverse function R[u) for u 
between and 1 imphes that u'{R{l)) > and u'{R) > for < i? < R{1). This 
immediately translates into the positivity conditions 

/(I) > 0, f{x) > for < a; < 1 (5.7) 

with /(I) = iff the model is critical in the sense of Section 3.2. Furthermore, we 
also obtain that p < 1, since u'{R) < 1 — GgRhj (^73|). Second, by the discussion of 
Section 4.1, we have r < 1 and the case r = 1 corresponds to a non-generic critical 
point. Note that this is precisely the range on which the equation ( |5.5| ) is well-defined. 

Indeed, since i7ik/16''){^^) ^Iforfc — )■ oo, the radius of convergence of il^{t) is 1, with 
'il'{t) ~ 1/(1 — t) for t — )■ 1, hence Kr{x, y) is a smooth function of (x, y) G [0, 1]^ except 
for r = 1 where it has a polar singularity at x = y = 1. The integral equation still holds 
for X = 1, which implies that /(y) — )■ for y — )■ 1 sufficiently fast. 



Conversely, given p and r both between and 1 and such that the equation (5.5) 
admits a solution / satisfying the positivity conditions ( |5.7| ), we may return to the 
original variables as follows. We first compute -R(l) via ( ^.6|) , then deduce u'{R), g 
and hi via (|5.4| ), and finally obtain u{R) as the primitive of u'{R) satisfying u{0) = 0. 
Following the steps of the above derivation backwards, we find that u{R) satisfies 



with the gkS given by (|2.3| ), so that we are indeed "solving" the 0{n) rigid loop model. 
The positivity conditions ensure that the inverse function R{u), thus the model, are 
well-defined. In particular, the generating function Fp°°^{n; g, hi,0) — Fp is directly 
expressed from / by 

'2p\ J^xPf{x)dx 
[f,'fix)dx) 

In conclusion, we have shown that generating functions for the 0{n) rigid loop 
model may be obtained by solving the linear integral equation (|5.5| ). In practice, the 
main difficulty is that the change of parameters {g,hi) — )■ (p, r) is rather intricate. 
Nevertheless, if we were able to compute the function / for arbitrary parameters n, p, r, 
then we could deduce F^°^{n]g^ /ii, 0) in a parametric form. Before further analyzing 
the solutions of (|0|), let us extend our formalism to the non-rigid case. 



23 



5.2. Extension to the non-rigid case 

We may repeat the above analysis in the case of arbitrary values of hi and /i2, 
replacing the fixed point condition ( |2.tj| ) by the more general one (ST^). The main 
complication is that an extra sum over a variable k' is involved. The equation ( [4 .81 ) 
allows to rewrite this sum as a contour integral. Furthermore by (|5.2|) we have 

F(.) = r 4mL. (5.9) 

Jo Vl-AzR 

Substituting these expressions into ( |2.2D , then into ( |5.1j ) , we may explicitly evaluate the 
sum over the variable k and derive a linear integral equation for u'{R). Again a more 
compact form is obtained after a suitable change of variables: we let x, y, f{x), p be as 
in ( p. 41 ) while we now define 

T = 4R{l){hi + 2h2). (5.10) 
After some work, we arrive at the same form for the linear integral equation 

f{x) + — Kr{x,y)f{y)dy=l-px (5.11) 
27r Jq 

but with a different kernel Kr{x, y), now given by 

f 2iC (1 - ryCy/^ 2\{l- rxA+(C))^/^ (1 - rxK.{Of/^ ) ^'''"^ 

where we may choose as integration contour the circle \C,\ = Co for any {1/ {tx)) < 
Co < 1/(7"?/), and where we introduced the rescaled functions 

hi + 2/i2 

Note that A_|_ and A_ are functions of the ratio /ii//i2, and therefore Kr^x.y) depends 
implicitly on this ratio. Despite the apparent complication in the expression for the 
kernel Kr{x,y), much of the discussion of Section 5.1 can be generalized to the non- 
rigid case. Using exactly the same arguments, we find that: 

- the consistency relation (|5.6|), 

- the positivity conditions ( ^.7|) , where again /(I) = iff the model is critical, 

- the inequality p < 1, 

- the expression ( |5.8D for the 0{n) loop model generating functions Fp°°P{n; g, hi, /i2) = 
p 

all still hold in the non-rigid case. The discussion of r, in view of that of Section 4.2 
and of the general expression ( |5.1C1| ), becomes slightly more involved. We nevertheless 
find that r is still between and 1, with t — 1 corresponding to a non-generic critical 
point. Furthermore, Kr{x,y) is a smooth function of {x, y) G [0, l]'^ except in the case 
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r = 1 where, remarkably, we obtain the same singular behavior as in the rigid case, 
provided that hi > 0. More precisely, when x and y tend to 1 (keeping (1 — a;)/(l — y) 
finite), we have 

Ki{x,y)r^-^. (5.14) 
1 — xy 

For the record, let us briefly explain how this property results from ( p. 121 ) by a saddle- 
point approximation. The contour integral in the latter equation is dominated, for 
X, y — )■ 1 by the vicinity of C = 1 and we set 

C = l + ieZ, x = l-eX, y = 1 - eY (5.15) 

with e — 7- 0. At leading order in e, we have 

A+iC) = l-e^* Z +■■■, with = -A;(1) = -A;(z*) = 1 (5.16) 

using again the "miraculous" involutivity of X+{z) around z*. We flnally obtain the 
estimate 

, , 1 p rfZ 1 11 1 , 

"^^^^'^^ ~ 7 y_ T(y-iZ)V2(x + iZ)3/2 = -exTY ~ " ^'-'^^ 

In the case hi = 0, this estimate must be doubled since A_ has then an equal, instead 
of negligible, contribution. We thus have Ki{x,y) ~ 2/(1 — xy) in this case. 

In summary, the generating functions Fp°°P(n; (7, /ii, /i2) for the 0{n) loop model 
on tetravalent maps may be expressed, via a change of parameters {g,hi,h2) — )■ 
(p, r, hi/h2), in terms of the solution of the linear integral equation ( |5.11|) . 

5.3. Discussion of the solution of i\5. 1 j| j and of its singular behavior 
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Fig. 9: Qualitative phase diagram of the 0{n) loop model in the (p, r) 
plane for flxed values of n (between and 2) and of /ii//i2- 
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As mentioned previously, Eq. (|5.11| ) is a Fredholm integral equation of the second 
kind for the function f{x). This equation depends on n, p, r and the ratio /ii//i2 
as follows: — n/(27r) is the so-called parameter of integral equation, p appears only in 
the right-hand side and r and hi/h2 determine the kernel ( |5.12| ) (in particular, for 
hi/h2 = oo, we recover the rigid case of (|5.5|) ). Here our terminology is borrowed from 
||iy|| . In this Section, we shall assume that n and hi/h2 are fixed quantities and will 
look at the dependence of f{x) on r and p, both varying a priori between and 1. We 
shall successively consider the case r < 1 and the case r = 1, since the latter is special 
and corresponds, as shown before, to non-generic criticality. Under certain reasonable 
assumptions, Fredholm theory implies the qualitative phase diagram of Fig. 9 in the 
(p, r) plane. This picture will be corroborated by the exact results of Section 6. 

The case r < 1: here Kr{x,y) is a smooth function of (x, y) varying in the domain 
[0,1]^, it is thus square-integrable. We may then apply Fredholm theory: assuming 
that — (27r)/n is not a characteristic value of the integral equation (i.e. there exists no 
non-zero function /(x) such that /(x) -|- (n/27r) Kr{x, y)f{y)dy = - we expect this 
assumption to be valid when n is between and 2), Eq. ( ^.11|) has a unique solution 
f{x) for each value of p. More precisely, by linearity, we have 

fix) = hix)-pUx) (5.18) 

where fi{x) (resp. fid{x)) is the solution of the linear integral equation obtained by 
changing its right-hand side into 1 (resp. x). Note that fi{x) and fid{x) do not depend 
on p but implicitly depend on r as well as n and hi/h2- These functions might be 
expressed, for instance, via Neumann series at least for n small enough. 

However, the solution f{x) does not necessarily satisfy the positivity conditions 
( p77|) ensuring that the model is well-defined. In view of (|5.18| ), we conjecture that these 
conditions amount to p < Pc, where Pc = /i(l)//id(l). In particular, when p = pc, 
/(I) = and the model is critical. We identify this critical point as generic in the sense 
of Section 3.2, since we expect /'(I) hence u"{R{l)) to be finite for any r < 1. As 
Pc depends on r, we obtain a generic critical line in the (p, r) plane, see Fig. 9. Note 
that, for any values of n and hi/h2, the generic critical line starts from (p = l,r = 0), 
corresponding to the critical point for pure quadrangulations without loops. Indeed, 
since Kq = 0, we have f{x) = 1 — px for r = 0, therefore the critical value of p is 1 (we 
also recover the known critical values -R(l) = 2 and g = 1/12). Moreover, for < ti < 2, 
we expect the line of generic critical points to connect continuously with the non-generic 
critical line r = 1 (to be discussed below): p^ should be positive for all r < 1, and have 
a positive limit p* as r — )■ 1. 

The case r = 1: now Ki{x, y) is no longer a smooth function of (x, y) in the domain 
[0, 1]^, but diverges as ( |5.14|) for x, y — )■ 1. We do not know whether a general theory 
applies to such kernels. Nevertheless, from |]2^, we expect that Eq.( [5.111 ) has, for all 
p, a unique continuous solution satisfying /(I) = 0. Thanks to this cancellation, the 
integral Ki{l, y)f{y)dy may still be well-defined. 
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Furthermore, we see that the solution f{x) cannot be regular for x — ?■ 1 (i.e. vanish 
as an integer power of x — 1) as otherwise, the integral in (|5.11| ) would contain singular 
terms (with logarithms) which are not present in the r.h.s, regular at a: = 1. This 
suggests to assume that the solution behaves as 



/(x)~C(l-x) 



(5.19) 



for some constant C and some positive non-integral exponent a (the condition a > 
ensures both that /(I) = and that the integral in ( 5.11 ) converges at x = 1). 
Remarkably, a cannot take arbitrary values, but is related to n via 



n 



2 sinyra 



as seen from the following argument. Subtracting to ( |5.11| ) its expression at x 
get 

f{x) + — / iK^ix,y)-K^{l,y))fiy)dy = p{l-x) 

'0 



2n 



(5.20) 



1, we 



(5.21) 



with Ki{x, y)—Ki{l, y) ~ —{l—x)y/{{l—y){l—xy)) when x and y tend to 1. For x — )■ 1, 
the dominant singular term in the integral arises from y's such that 1 — y = 0(l — x). 
More precisely, writing y = l — s(l — x),we get at leading order 



(i^i(x,y)-Ki(l,y))/(2/) dy 



-C(l -x) 



sing. 



1 + S 



-ds 



(5.22) 



with the right-hand side integral evaluated as tt/ sinTra. As such, this equation holds 
only for o: < 1 so that the integral in the r.h.s is convergent. It remains valid for a > 1 
provided the integral is analytically continued, with again the value tt/ sinrra (note that 
the integral in the l.h.s also contains regular terms that, for a > 1, dominate the leading 
singular term). Combining with (|5.19D back into ( ^.21|) , whose r.h.s contains no singular 
term, we deduce ( |5.20| ). 

Note that, for a given value of n between and 2, the possible a > satisfying 
( ^.20|) are of the form a = 1/2 ± 6 + m, with 7fb = arccos(n/2) and m a non-negative 
integer. So far we have only discussed the dominant exponent but, by linearity, all 
exponents appearing in the expansion of /(x) at x = 1 should be of this form. We 
therefore expect /(x) to be of the general form 



fix) 



X 



,1/2-b 



$_(x)-f(l-x)i/'+''$+(x) 



(5.23) 



where $_(x) and $+(x) are functions with a regular expansion at x 
only integral powers of x — 1 



1 (i.e. contain 

By considering the expansion of (|5.11| ) at x = 1 and 



splitting it into singular and regular parts (i.e. separating terms with non-integral and 
integral exponents), we find that $_(x) and $+(x) should satisfy 



{1 - xy/^^^^±{x) + 



cos(7r6) 



TT 



Ki{x,y){l-y)'/^^'^±{y)dy 



(5.24) 



sing. 
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and 



1 — px 



reg. 

(5.25) 

Now we expect that $-(1) is not zero generically so that the leading singularity at 
a; = 1 is of the form (1 — x)^/'^~^. Furthermore, $-(1) must be non-negative in order 
to satisfy the positivity conditions ( ^.71) . This situation should hold for any p < p* 
and corresponds to the dense phase of the 0{n) loop model, see Fig. 9. At p = p* , we 
expect $-(1) = and $+(1) > so that, at this special point, the leading singularity 
becomes of the form (1 — x)^/^"'"'', corresponding now to the dilute 0{n) loop model. 
For p > p*, we expect $-(1) < so that the model is ill-defined. That the critical 
condition /(I) = for r < 1 coincides when r — ?■ 1 with the condition $-(1) = is 
rather natural and makes us believe that the transition from the line of generic critical 
points to that of non-generic ones should be continuous. Note that having an effective 
value of a larger than 1 for the leading singularity simply corresponds to a situation 
where both $-(1) and $+(1) would vanish. This may occur only for a particular class 
of right-hand sides in the integral equation and it is not expected in the present case. 

Finally, recalling Section 3.3, we see that the behavior ( |5.19|) for /(x) = u'{R) 
corresponds precisely to the expansion ( |3.12| ) for u{R). In particular, we identify 

3 

a = CK + - (5.26) 
2 



and the relation (|5.20|) is nothing but ( [4.4|) . 



6. Exact phase diagram of the rigid loop model 

In this Section, we concentrate on the rigid loop model and explain how to derive 
its exact phase diagram. 

6.1. Equations for the resolvent 

The starting point is a linear integral equation, not for the function u'[R) as in 
Section 5.1, but for the resolvent W{^) defined as in ( |3.19| ). Those two quantities are 
related via the general formula ( p.20| ) . Moreover, substituting the fixed point condition 
I) for the rigid loop model into (|3.4| ), we have for R < R{1) 



u{R) =R- 3gR^ - ^ (^^) hfR^F, 

(6.1) 



R-,,R^-^-l'-^[ ^ ' 

2 J 2m \^l-4hlRe 
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where the contour of integration is, say, the circle of radius 7 = 2a/-R(1) (note that 
4/i^i?7^ is always smaller than 1 by Eq. ( [4.2|) ). Upon differentiating with respect to R, 
and substituting into ( p.20| ), we obtain 



-3/2 



2^ / 2i7r hli'"^ - 1/^2 \^ Y 1 - 4i2(l)/ife2 

where 



(6.2) 

s{i) = \[i-ge- i{i - 2^7i^(l) - gewi-m^)/e) (6.3) 



corresponds to the first two terms in the r.h.s of ( |6.1| ). 

Eq. (^]^) is a linear integral equation for the resolvent which rephrases that of 
Section 5.1 for /(x) = u'{R). It implies a simpler functional equation for W{^) as 
follows: for ^ G [—7, 7], we have 

S{i + iO) + S{i - iO) = e - ge (6.4) 

so that, from ( |6.2| ), we may write 

The latter integral may be evaluated by the residue theorem: the integrand has poles 
at ^' = ±(/ii^)~^, each with residue {hi^)~^W{{hi^)~^) /2, and a pole at ^' = cxd with 
residue —Fq = — 1. Hence the resolvent satisfies the functional equation 

Til Til / 1 \ 

+ iO) + _ iO) = ^ - ^ - - -^W { — \ , e G [-7, 7]- (6.6) 



Note the similarity with Eq.( p..2D for the different model of Section 1.2. Equation 
can also be obtained as a consequence of loop equations in the matrix model formulation 
of our 0{n) model [|l], Eq.(V-22)]. 

Solving ( |6.2p boils down to finding a solution of ( |0| ) which is bounded, odd in ^ 
and such that VF(^) ~ 1/^ as ^ — 00. As we shall see, these requirements fix W{^) 
completely. By linearity, we may write 

W^(0 = W^part(0 + W^hon.(0 , (6-7) 

where VFpart(0 is the easy particular solution of ( |6.6|) 



^part(0 = J ^ + T^T^, (6.^ 

4 — (2 + n)4 
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and where Whom(0 is now an odd solution of the homogeneous equation 

W^hom(e + iO) + M^hom(e - iO) + /^^hom (^^^ = . (6.9) 

The condition that W{^) ~ 1/^ for ^ — )• oo, and that W is bounded amounts to 
demanding that 



77,2 4 — 77,2 2 + n " 

/ 2o 2 2 \ ^^-^^^ 

We shall give in Section 6.3 the general expression for Wh_om(0 when 4hiR(l) < 1, which 
corresponds to a subcritical or generic critical situation. Something special happens on 
the non-generic critical line 4hiR{l) = 1 because then, the cut [—7,7] collides with its 
image under ^ 1— )■ l/{hi^) (see Fig. 10 below). As we shall now see, V1^hom(0 ^ simple 
expression along this line. 

6.2. Non-generic critical line 

As seen in Section 4.1, the non-generic critical line is characterized by r = 
4hiR{l) = 1 (or equivalently -R(l) = Rc)- This implies that the extremity of the 
cut 7 = 2^R{1) = l/^/fH is a fixed point of ^ 1— )■ l/(/ii^). It is then possible to guess 
the general solution of the homogeneous equation ( |0| ) which is odd in ^, namely: 



(6.11) 

where B{^) is some arbitrary analytic function and where, again, nb = arccos(n/2). 
That this form satisfies (|6.9| ) can be checked directly and one can prove that it describes 
all the solutions. This in turn leads to a spectral density supported on [—7,7]: 





(6.12) 

The requirement that Whom is holomorphic in C \ [—7, 7] imposes that B is an en- 
tire function. To satisfy ( p.lO|) , B{^) must be a polynomial of degree 3, whose four 
coefficients are determined from 

^^^)f7T^y = r^^'-r^^+^^"' + ^^^"') as ^-^00. (6.13) 
Vs+7/ 4 — 4 — 77,^ 2-|-n 

Note that we could a priori imagine terms of order and in this expansion since 
they would be canceled by parity in the expansion of VFhom(0 at large ^. However, 
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such terms would create poles for p(^) at ^ = 0, which are not allowed because p should 
be integrable. The condition ( |6.13|) amounts to five equations: four of them fix the 
coefficients of ^^'^ the last one yields some additional relation between g and hi 
which is nothing but the equation for the non-generic critical line. We find explicitly 



^(0 = ( + 267^2 ^ 262^2^ + ^(6 + 2}f')A - —^{^ + 267) (6.14) 



while the equation for the non-generic critical line reads 

3 /, 2-n 



In order for the expression ( |6.11| ) to be consistent with the positivity of u'{R), the 
spectral density ( |6.12|) must be positive on ] — 7,7[. In particular, expanding p(^) for 
^ — > ±7, we have 



,«) = !i^^M('4'(B(-,)-,B'(-,))(l-S;) 



+4-''(S(7)+7S'(7))(l-^J 1+0^1 



(6.16) 



Demanding the positivity of p{^) requires i?(— 7) — 7i?'(— 7) > 0, which yields the 
condition 

3hi 

' ^ 2(6^-26 + 3) ■ ^'-''^ 

The non-generic critical line ( |6.15|) cannot extend outside of the region defined by ( |6.17|) 
and therefore ends when it hits the boundary of this region, i.e. at a point {g*, h^) with 

36^(2-6)2 6^(2-6)2 
^ 2(2-n)(62_26 + 3)2 ' ^ (2 - n)(62 - 26 + 3) ' ^ ^ 

As shown in Section 6.4, the critical line becomes generic beyond that point. When 
{g,hi) = {g*,hl), we may check that the coefficient of the first subleading term is 
positive, i.e. -6(7) -|-7-B'(7) > 0. Comparing ( |6.16| ) with (|3.23| ), we therefore read the 



value of the exponent a, namely a = 2 — b for g < g* and a = 2 + b8itg = g* 
6.3. General expression of the resolvent 



Eq. ( |6.6|) can be recast more elegantly in terms of the differential form w(^) = 
VF"hom(0 di-: namely: 

G [-7, 7], uj{i + iO) + uj{i - iO) - nu{s{i)) = , (6.19) 
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oo 




Fig. 10: The cut of W{^), or equivalently of WhomiO (solid line) and its 
image (dashed line) under the involution ^ l/{hi^). 
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Fig. 11: Construction of the elliptic parametrization v{^) of Eq. (6.20), 
which depends on the path followed from — (/ii7)~^ to ^. The upper and 
lower half-planes map respectively to the left and right rectangles (whose 
union is denoted by V). 



where s is the involution ^ l/{hi^). Note that (s(^))^ = 1/A+(1/^^), where A_|_ is 
the involution of Section 2.3, specialized to the rigid case. We underline the similarity 
of Eq. (|6.19| ) with that relevant in the 0{n) model discussed in Section 1.2, where loops 
visit only vertices of degree 3: in this case, we had a different involution ^ i— )■ — ^. 
So, the techniques already developed for the 0{n) model where loops visit only vertices 
of degree 3 [jll|,^ can be applied to Eq. ( |6.19|) with few modifications. 

The general strategy is to introduce a parametrization which opens the cut [—7, 7] 
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as well as its image ] — oo, — (/ii7) ^] U [(/ii7) ^, +oo[ under the involution (see Fig. 10), 
for instance: ^ 

v{^) = V [ "^"^ . (6.20) 

The new variable v can be expressed in terms of Jacobi elliptic functions and is a 
multivalued function of ^, depending on the path followed in the complex plane from 
the origin (/ii7)~^ to ^. Conversely, one may view ^ as a function of v e C. We choose 
the constant V by demanding that v{{hi'-y)~^) = —1/2 when following a path with small 
positive imaginary part from — (/ii7)~^ to oo and then back to {hi^)~^. This leads to 

where K is the complete elliptic integral. We then denote T — f (— 7) (by the most 
straightforward path), and the fact that the square root discontinuity is included in the 
real axis implies that T = i\T\. The function ^{v) is easily seen to have the following 
properties: 

from which one can deduce also the properties of its derivative $,'{v). For bookkeeping, 
we mention the expansion when v — > Voo = —1/4: 

^v) = + E^{v- v^) + 0{{v- v^f) (6.23) 

V - Voo 

with 

•T/ -1 

= iV 



SiS_i = i(7^ + (/zi7)-^) , 

where we denote K'{t) = i^(vT^-^). The image of the points of C \ ([— 7,7]U] — 
oOj —{hi'~f)~^] U [(/ii7)~"'^, +00Q reached by a path which does not cross the seg- 
ment [7, (/ii7)~^] (for instance a straight path) is the domain V — {v : Re[f] G 
] - 1/2, l/2[ ,lm[v] e]0, |T|[}, as shown in Figs. 10 and 11. Let us define 

rz;{v) = W^hom(e(^))^'(^) (6.25) 

which is an analytic function on V. Since Whom(0 discontinuity when ^ e 

] — 00, —{hi'y)~^] U [(/ii7)~^, oo[, w{v) takes opposite values when v ■<->■ — v along the 
segment [—1/2,1/2]. Thus, w can be extended to an analytic function defined on 
V2 = V U (-V) by setting 

e -V, w{v) = -w{-v) . (6.26) 
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Likewise, the absence of discontinuity along ^ G [7, (/ii7) ^] allows to extend w as an 
analytic function defined on the strip Vstrip = Umez (^2 + ^) by setting 



G (V2 + m), tA7(f ) = ti7(f — m) . 



(6.27) 



The top boundary of this strip, Vcut = \y '■ Ini[f] = maps to points ^ G [—7,7]. 
Eventually, looking at Eq. (|6.22D and using Eq. (I^H), Eq. ( |6l9D turns into: 



Vi;GVcut, + - 2T) - - T) = . 



(6.28) 



This allows to extend recursively and without ambiguity w as an analytic function on 
^ — UmezC'^strip + 2mT). For instance, for all v G Vstrip + 2T, w{v) can be defined as 
nw(y — T) — ro(f — 2T), noticing that v — T and v — 2T belong to Vstrip- Since (|6.26| ), 
( |6.27|) and (|6.28|) are linear analytic relations between analytic functions, they are now 
valid for any value f G C. 

To summarize, we reduced the problem to that of finding an analytic function 
V I—)- w{v), which is odd, 1-periodic, and satisfies: 



w{v — 2T) — nw{v — T) + vj{v) 
Eq. (|6.10| ) demands that ro(f ) behaves as: 







2gEU 



vj{v) ■■ 



1 



2El, 



A — 'n? [v — Vr 



l-f(7^ + (M)-^) 



{v — VooY 2 + nv — Vc 



(6.29) 



-0(1) 



when V ^ Vr 



-1/4, and: 



zu{v) ■■ 



"4 



1 



n 



^2 ^ 



A — n^ (f — fo) A — n^ 



l-|(7^ + (/^i7)-^) 



1 



n 



1 



- (v — vq)^ 2 + n v — vq 



(6.30) 



■fO(l) 



when V ^ vq = T — 1/A. We now assume n i {0, 2}, such that e^^^^ ^ 1. Defining 



(6.31) 



g±i7rb g=pi7rb 



(6.32) 



we have w{y) = w^iy) + w-{v)^ while, from the 1-periodicity of w and from (|6.29|) , 
vj± must satisfy: 



(6.33) 



Such functions are generalizations of elliptic functions, and they can be constructed by 
taking appropriate ratios of the Jacobi theta function of nome q = e^^^. Let us just 
state the existence of a unique analytic function (h, which satisfies 



a{v + 1) = Cb{v), Cb{v + T)= e'^'Cbiv) 



(6.34) 
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and has a unique pole when v = mod Z © TZ , which is simple, and is such that 
Cbiv) ~ 1/f when f — )■ 0. Its construction and main properties are listed in Appendix 
B . One may generate functions satisfying ( |6.33p with poles of higher degree at f = by 
considering the derivatives of or of C-b{v) = —(bi—v), and put this pole at any given 
point w by shifting the argument v to v — w. Since holomorphic functions satisfying 
( |6.33| ) must vanish identically, one may determine zu by matching the divergent behavior 
at its poles with a linear combination of the previous functions. This leads eventually 
to: 



zu{v) 



1 



2 + n 



V {C,{v - 1/4) + Cb{v + 1/4) - ai-v + 1/4) - a{~v - 1/4)} (6.35) 



where V is the differential operator: 
V 



24(2 -n) " 2(2- 
Then, the spectral density is given by 



n] 



-l + f(7^ + (M)-^) 



(6.36) 



2i7r 
dw 



(I^hom(^(2T - V)) - Wi^omiCiv))) dv 
1\ . / 1~ 



-n 



2 + n 



VUb w 



1 



+ Cb( w + - ]+Cb(-w + - ]+Cb(-w- 



1 

4 

(6.37) 

where we have set w = T — f so that w G [0,1/2] corresponds to ^ G [—7,7]. The 
value of 7 = 2a/-R(1) is determined a posteriori as a function of g and hi hj requiring 
from (|3.21| ) that p(^) oc ± 7 when ^ — )■ =1=7, which is equivalent to demanding that 
p{C{v))C'{v) = 0{w^) when w ^ 0. 

6.4- Generic critical line 

In the solution above, the generic critical line is the relation between g and hi 
obtained by demanding that p(^) oc (^ ± 7)^/^ when ^ — )■ =F7 (see Eq. ( |3.23| )). In other 
words, in the Taylor expansion of p{^{v))^'{v), as given by ( |6.37| ), when w — )■ 0, the 
generic critical line is characterized by the vanishing of the terms of order 1 and w'^. 
We may write these two conditions in a parametric way, with parameter r = /ii7^ = 
4/iii?(l), as 



Zo + 
Z2 + 



where 



^2 ^ 


(- 
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2(2 -n) 


3hi 


^2 _^ 


(- 
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2(2 - n) 


3/ii 


Z2j = 







(r + r-i; 
(r + r-i; 



Z2 + 



^4 + 



^U9 



n] 



24(2 
24(2 - n) 







(6.38) 



(6.39) 
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Note that depends implicitly on r via C& which depends on T, itself related to r via: 

,^?3(0|4T)y' ^3(01^) 

The solution is: 

1 As 



^(^.^}^\ . (6.40) 
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2(2 - n)r [K'(r)]2 Ag + 4{t^ + 1) [K'(r)]2 A4 

6 Ag A4 

2^ (As + 4(r2 + 1) [K'(r)]2 A^)^ ' 



(6.41) 



where K'(t) = K(\^l — r^) and the A's are functions of r defined via: 

A4 = T\ZoZ^ - Z|), Ae = -T6(ZoZ6 - Z2Z4), Ag = T^iZ^Z^ - Zj) . (6.42) 

Those equations become simpler in the neighborhood of the special points r -^■ 
(i.e. /ii — )■ 0), or r — )■ 1 (i.e. hi — )■ hi, the tip of the non-generic critical line). When 
hi —7- 0, it is convenient to use the variable q = e"^^ — ?■ for asymptotics, and we find: 

2 2 (6 43) 

/^i = y-(2+^)?^+(7 + 2n+^) q'' + 0{q'^) 

As expected, the generic critical line meets the critical point of quadrangulations at 
{g = 1/12, hi = 0), and near this point it behaves as: 



When r —7- 1, it is convenient to use the variable q' = e"''^/^ — > for asymptotics. 
The computation confirms that the generic critical line ends at the tip {g*, hi) of the 
non-generic critical line found in Eq. (|6.18|) , and we find near this point that: 

g-g*\ 2(1-26) f hi - hl\ (2-6)2 f hi - hi 



g* J 2 + \ hi J 2 + \ hi 

64(1 + b) / (2 - 6)2(1 - 6)2(3 - 26 + 62) \ f hi - hl^ 



^ (l-6)(2 + 62) V 46(1 + 6)2(2 + 62) J K hi 

(6.45) 

The first two terms coincide exactly with Eq. ( |6.15|) : when passing from the generic to 
the non-generic critical line, both the slope and the curvature remain continuous. For 
n = 1 (i.e. 6 = 1/3), we find a leading discontinuity in the third derivative of g with 
respect to hi , as expected for the Ising model ||2^ . 

36 



h 



1 0.30 



0.25 



0.20 



0.15 



0.10 



0.05 



non-generic 
^ critical 



b=0.3 



dilute^ 



subcritical 



generic 
critical 



0.00 



0.00 



0.02 



0.04 



0.06 



0.08 



1/12 



0.10 



9 



Fig. 12: The exact phase diagram of the 0{n) rigid loop model in the 
{g,hi) plane. It is shown here for h — 0.3 (with n = 2cos7r6) but it is 
qualitatively the same for any value of h between and 1/2. The critical line 
separates the region where the model is subcritical from the region where 
it is ill-defined. The type of criticality changes along the line: generic for 
g > g*, non-generic for g < g* with an exponent a = 2 — b (dense model), 
non-generic at g = g* with an exponent a = 2 + b (dilute model). The line 
of non-generic critical points is an arc of parabola which we extended in 
dashed line for clarity. 



6.5. Phase diagram 

The results above are best summarized in the phase diagram of Fig. 12. We have 
found a line of non-generic critical points given by the arc of parabola ( |6.15| ), which 
hnks the point {g = 0, /ii = 26^/(2 - n)) to the point {g*, hi) of Eq. ( |6T8|) . Along this 
line, the exponent a takes the value 2 — 6. At the terminating point {g*,hl), a takes 
instead the value 2 + b. We then found a line of generic critical points with a more 
complicated parametrization ( |6.41| ). When r decreases from 1 to 0, this line links the 
point {g*, hi) to the point {g = 1/12, /ii = 0) describing pure quadrangulations. As just 
mentioned, the non-generic and generic critical lines connect with a continuous slope. 
Their concatenation forms the line hi = hc{n]g) of Section 4.1 as the model cannot be 
well-defined above this line. 

Let us finally note that, in the limit n — )■ 0, then g* — )• 1/12 and hi — )■ 1/8. The 
non-generic critical line tends to the arc of parabola 

g=^{hi-4hl), ^</ii<^, (6.46) 
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Fig. 13: The exact phase diagram of the rigid 0{n) model in the limit 
n 0. 

while the generic critical line becomes the vertical segment parametrized hy g = 1/12, 
< /ii < 1/8 (see Fig. 13). This may be understood as follows: the small n expansion of 
describes quadrangulations equipped with a fixed finite number of rigid loops. Us- 
ing exact enumeration results for quadrangulations with multiple boundaries, it can be 
seen that the contribution to -Fp°°P from loops of large length 2k behaves as {AhiRq)'^^ ^ 
where Rq is the generating function -R(l) for pure quadrangulations, solution of 



(6.47) 



Having a finite contribution from large loops requires that AhiRq < 1. Note that Rq 
ranges from 1 to 2 when g ranges from to 1/12. Criticality may be obtained in two 
manners: either we set Ah^Rq = 1, so that the contribution from large loops decays 
sub-exponentially. Note then the equivalence between Eqs. (|6.46|) and ( |6.47| ). Or we 
may set g = 1/12, so that the pure quadrangulations are themselves critical. This last 
situation requires hi < 1/8, since Rq = 2 in this case. 



7. Variants of the 0{n) loop model 

In this Section, we briefly discuss other versions of the 0{n) loop model on quadran- 
gulations, including models with non-symmetric local weights and models with restricted 
loop lengths. More precisely, we concentrate on non-generic critical points and discuss 
how the relation ( [4.4|) is modified in these cases. We finally extend our results to maps 
whose faces have arbitrary (but bounded) even degrees. 
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7.1. N on- symmetric models 



At this stage, it should be clear to the reader that the relation ( |4.4| ) between the 
loop weight n and the exponent a only depends on a few properties of the ring transfer 
matrix M{z). Denoting by \+{z) the largest eigenvalue of M^(z), we used the estimate 



-k' 



X^(z)'^ to obtain eventually the relation 



n 



2(;U*)"sin7r(a-3/2), ^* = -A' 



(7.1) 



where z* is the fixed point of the mapping z i-> X^{z). In the symmetric case discussed 
so far, this mapping is an involution in the vicinity of z* , so that X'_^_{z*) = —1. Let us 
now consider a slightly modified, non-symmetric version of our model defined as follows: 
the squares of type (c) in Fig.l come in two species, those whose two edges not crossed 
by the loop belong to the inner contour, and those where these two edges belong to 
the outer contour. We may as well view these two species as corresponding to outward, 
resp. inward turns of the loop at hand. For instance, the loop in Fig. 3 makes 8 inward 
and 6 outward turns. Assigning now different weights, say /i2,out and /i2,in respectively 
to these squares, the transfer matrix M{z) is replaced by 



M(z) 



I hiz-y^ 



1 



\ 



1 — ^2, out 
^2, in 



1 — h2,outZ~ 





(7.2) 



The relation (p.7| ) becomes 



2 9 

(/i2,outA + /l2,in^ - A^) - h-^Xz = 



(7.3) 



whose largest solution Xj^{z) now has a fixed point at z* = hi + /i2,out + ^2, in- Since 
A and z do not play symmetric roles in (|7.3|) , the mapping z X+{z) is no longer an 
involution and we now have a non-trivial value //* = {hi + 2/i2,out)/(^i +2/i2,in), leading 
to the new relation 

Note in particular that the range of n allowing for non- generic criticality is modified. 
It is instructive to compare this result to that obtained on a regular square lattice. In 
this context the parametrization n = 2 sin n{a — 3/2) = 2 cos Tcb naturally appears in the 
Coulomb gas approach to the model [23|. Many critical exponents of the model have 
simple (typically polynomial) expressions in terms of a (or b) . On such a regular lattice, 
we also have a well-defined notion of exterior and interior of a loop and we may give 
different weights to outward and inward turns. On a regular lattice, there are however 
4 more inward than outward turns, so that the symmetry of the model may be restored 
at the price of a rescaling n — > n(/i2,in/^2,out)^ of the weight per loop. This in turn 
changes the relation (|4.4|) into n = 2sin7r(a 



3/2)(/l2,out//l2,in)^. 
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Returning to the 0{n) model on quadrangulations, the case hi = is, as before, 
special since the largest eigenvalue is degenerate in this case, hence (|7.4|) is replaced by 

n = sin7r(a-3/2) f^^y . (7.5) 

7.2. Loops with restricted lengths 

Returning to the symmetric case, we may impose some restriction on the lengths 
of the loop by demanding for instance that they be multiples of a fixed integer A^. 
Such a restriction may occur for instance when the loop model is inherited from some 
underlying edge coloring problem. Now the length of a loop whose outer and inner 
contours have length 2k and 2k' is simply [k + k'). The consistency relation ( |2.2| ) has 
to be modified to account for the new constraint, and we are naturally led to consider 
now the quantity 

N-l 

fc'>o j=0 

fc + fc'=0 mod JV 

(7.6) 

to be estimated as before in the vicinity of z* and for large k. The above sum behaves 
as (1/A^) (A_|_(z*))'^ with now a prefactor provided \\j^{z*u~^)\ < A_|_(z*) when 
j = l,---,A^ — 1, which holds for /i2 > (again we also suppose that /ii > to avoid 
that \-{z*) = \j^{z*)). The correcting factor (1/A^) trivially results into a change of 
the relation ( |4.4| ) into 

n = 2A^sin7r(a-3/2) (7.7) 

(for /ii, /i2 > 0). In particular, imposing an even size for the loops takes the 0{n) loop 
model in the universality class of the 0{n/2) loop model without the parity constraint. 
This fact was already recognized in Ref. ||2^ in the slightly different context of loops 
living on triangles. 

In the rigid case /i2 = 0, loops are automatically of even length by construction. 
We have \j^{u~^ z*) = u^^ X^{z*) and the above sum behaves as A_(_(z*)'^, as before, 
provided /c is a multiple of for odd (respectively a multiple of for even) while 
it vanishes otherwise. For /i2 = 0, the relation ( |4.4| ) is therefore unchanged. Finally, for 
hi = 0, we get n — N smn{a — 3/2) instead. 

7.3. Faces with arbitrary even degrees 

Our results are easily extended to the case of maps whose faces have arbitrary even 
degrees, provided these degrees remain bounded, say by 2M + 2 (M > 1). Faces of 
degree 2m not visited by a loop receive a non-negative weight g^'^^ (1 < m < M + 1) 
and those visited by a loop receive a weight /),("^i'"*2) [jxii,m2 > 0) if the face has mi 
(resp. m2) incident edges belonging to the outer (resp. inner) contour of the loop at 
hand. Since the total degree of such a face is mi + m2 + 2, we will implicitly assume 
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in the following that h^^^'^^^ is non zero only if mi + m2 < 2M and mi and m2 have 
the same parity. With these new weights, the 0{n) loop model is now described by the 
fixed point condition 

M+l 

9k =Y^ 9^'"^Sk,m + nJ2 Ak,k'Fk' (7.8) 

m=l k'>0 

where Ak^k' is the generating function for (rooted) rings (now made of faces with arbi- 
trary even degrees) with sides of lengths 2k and 2k' . As before, this generating function 
is best encoded in the quantity 



k'>0 



k'Z 



-k' 



trM{z 



,2k 



(7.9) 



involving a new transfer matrix M{z) of size 2M x 2M now given by: 



M{z) 



1-^ /l{0.'"2)2-T' 



■m2) y 5^ 



— i ^ 

1-22 

1 

1 



\^f^i2M,m2)z-^ 






















(7.10) 



0/ 



Here the sum in 'Y^h^^' 



\ ^9 

m2) y Y- 



runs over m2 > (and in practice over values of m2 



ranging from to 2M — i and having the same parity as i). As shown in Appendix A, 
the eigenvalues A of M'^{z) are the solutions of the characteristic equation 



{Xz 



, , . . 2M — rriT 2M — m.o 



m-^ ,m'2 even 



m, -[^ , m2 ^ 
in A ,T7-io odd 



mi,m2) 



(7.11) 

Repeating the analysis of Section 4, we again find at a non generic critical point the 
consistency relation between n and the exponent a characterizing the large k asymptotics 
oiFk-. 

n = 2(;U*)"sin7r(a-3/2), ^u* = -A+(2*) (7.12) 

in terms of the largest eigenvalue A-|-(2;) of M'^{z) and its fixed point z* . Note that this 
relation holds when the largest eigenvalue is not degenerate. In a symmetric model, 
we must set /i(™'i''^2) _ }^im2,mi) that A and z play symmetric roles in ( [7.11| ). This 
implies as before that the mapping z — > X+{z) is an involution in the vicinity of z* and 
that ^* = 1 so that the simple relation (|4.4|) is recovered. 

This generic relation (|7T|) is modified whenever /i("^i'"^2) _ q fgj. q]\ yg^[_ 
ues of mi and m2. In this case, the r.h.s of Eq. ( |7.11| ) vanishes and the largest 
eigenvalue is degenerate, resulting in the suppression of the factor 2 in (|7.12|) , namely 
n = (//*)" sin 7r(a — 3/2) with again //* = 1 in the symmetric case. 
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8. Conclusion 



In this paper, we have shown how to relate a number of 0{n) loop models to models 
of bipartite maps. More precisely, we have shown that the gasket of an 0{n) loop model 
configuration is distributed according to a Boltzmann ensemble of bipartite maps with 
appropriate degree dependent face weights {gk)k>i- Those weights are determined by a 
fixed point condition inherited from a bijective decomposition of the 0{n) configurations 
along the contours of their loops. In particular, the non-generic (dense and dilute) 
critical points of the 0{n) loop models correspond to ensembles of bipartite maps with 
large faces belonging to the class considered in [l^] with a distribution characterized by 
some exponent a between 3/2 and 5/2, related to n generically via 

?i = 2sin7r(a-3/2) . (8.1) 

Technically, this formula is one of a number of simple consistency relations dictated by 
the fixed point condition at a non-generic critical point. Their derivation involves only a 
few properties of a simple transfer matrix M{z) describing the sequence of faces visited 
by a loop (the ring). For instance, the constant z* characterizing the exponential decay 

of the face weights {gk ~ (^*)^) or the exponential growth of Fj^'^^ (-^fc°°^ ~ (1/^*)'^) is 
simply obtained as the solution of the equation z* = A+(z*), where X+{z) denotes the 
largest eigenvalue of M^(z). 

Noticeably, the same scheme appears to work also for 0{n) loop models where the 
loops visit only trivalent vertices or, equivalently on the dual, where the ring is made 
of a sequence of triangles. The ring transfer matrix reduces in this case to a scalar 
M{z) = iX+{z)), with 

X+{z) = h/{l-h/z) (8.2) 

(note that we do not square M as contours are not required to have even lengths). Here 
again, the exponential growth factor {1/z*) = {2h)~^ for Fp°°P in ( |1.4| ) is the solution 

of z* = X+{z*), while equation ( |8.1| ) still holds. 

A corollary of our reformulation is that the metric properties of the gasket at a 
non-generic point may be obtained from those of ensembles of bipartite maps with large 
faces. In Jl^, it was shown that these maps have a fractal dimension 2a — 1 and one 



may hope to be able to extract the average gasket profile from known expressions for 
discrete distance dependent two-point functions in bipartite maps. 

Acknowledgments: The work of G.B. is partly supported by the ANR project 
GranMa "Grandes Matrices Aleatoires" ANR-08-BLAN-0311-01. 



Appendix A. Eigenvalues of the transfer matrix for arbitrary even degrees 

Here we consider the 0{n) loop model on maps with arbitrary even face degrees, 
as defined in Section 7.3. Let us introduce the quantity 

k>l k'>0 
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where Ak^k> enumerates configurations of rooted rings witli outer and inner sides of 
lengths 2k and 2k', with the face weights of Section 7.3. Due to the factor l/(2/c), S 
enumerates unrooted ring configurations of arbitrary side lengths, with a weight 
(resp. z~^^'^) per edge of the outer (resp. inner) contour. By a direct calculation, we 
have 



fc>i 

-^trlog(l -m;"^M^(2)) 

--logdet(l -«;-^M2(;2)) 
2 



(A.2) 



where M{z) is the transfer matrix defined in (|7.10|) . On the other hand, if we denote 
by i the length of a ring, i.e. its number of faces or equivalently the length of the 
underlying loop, the configurations counted by S are simply cyclic sequences of length 
i made of the various squares at hand. We may therefore write 



i>l mi,m2>0 



Z 2 



7712 

w~~ z~~ 



mi ,m2 >0 



e>l m2>0 

(A.3) 

Here the first two terms differ only by the factor (— 1)"^^, so that their half-sum selects 
configuration where the total length of the outer contour is even (and so is that of the 
inner contour since mi and m2 have the same parity whenever /i("^i'"^2) The third 

term is the — )■ limit of the first two and is subtracted to account for the fact that 
there is no /c = term in S. Summing over £, this leads to 

(i- J2h^"'^''^^H-'^)'^'w-^z-^^\ 

mi,m2>0 mi,m2>0 



^ = -llog 
2 ^ 



1 m-^ rn,2 

'^^>w~~ z~~ 



V m2>0 / 

(A.4) 

from which we readily extract the value of det(l — w ^M'^{z)) by comparison with (|A.2| ). 
We end up with the final expression 

P{X,z) 



det{X -M^{z)) 



Q{z) 



(A.5) 



where P(A, z) is the polynomial: 

P{X,z) = ({Xz)^' - 



/^(mi,m2)^ 2 



2M— mi 2Af — mo 



mi ,m2 >0 

mi,m2>0 



(A.6) 
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and Q{z) the polynomial: 



m2>0 



Writing P{X, z) — leads precisely to Eq. 

Appendix B. Properties of the function (i, 

The (first) Jacobi theta function •di{v\T) is defined by 



mez 



The series is absolutely convergent when Im[T] > 0, and the following properties can 
be easily proved: 

^i{v + l\T) = -^i{v\T), ^?l(^; + T|T) = -e'^^-2^^"i?i(^;|T) . (B.2) 

The point v = turns out to be the unique zero (modulo Z © TZ) of 'di{v\T). This 
function has the following modular property: 

^i{v\T) = ^^e-'^^i(^\^] (B.3) 



v/=iT \T' T 

which is useful to relate the limit \T\ — )■ to the limit \T'\ — t- oo, where T' = —1/T. 
Notice that when \T'\ — > +oo, we have q' = e^^^ — )■ 0, and therefore the series ( [B.l|) 
gives in a straightforward way the asymptotics of 'i?i(f |T'). Other Jacobi theta functions 
that appear in the text are: 

'&2{v\T) = ^ gi7r(m-l/2)2T+i7r(2m-l)^ 

^3(^^|T) = e'™'^+2^-^ (B_4) 



meZ 



ITT m'^T -{-2177 mv 



We also introduce the Weierstrafi elliptic function: 

(Z,m)6Z2\{(0,0)} ^ ^ > V ^ > 
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The function p is even, 1 and T periodic, and has the following properties: 

Ratios of -i?! may be used to construct functions having prescribed poles and zeroes, 
and taking a constant or linear phase when v ^ v + T. For instance: 

^^(^^ - Mv\T) M-b/m ^ ^ ^ 

is the unique function which is 1-periodic, takes a phase e^^^ when v ^ v + T, has only 
a simple pole at t; = OmodZ © TZ, and is such that Cb{v) ~ l/v when — > 0. The 
value of the phase under translation by T implies that has a unique zero (modulo 
Z © TZ), located at f = 6/2. We may also find with another representation: 

(^^(v) = J2 e-''^''^ TT cotan 7t(v + mT) , (B.8) 
mez 

where Xlmez ■ • • has to be understood as limM->oo (Zlm=-M ' ' ') • 
Let us introduce coefficients Cq and C\ such that 

Cb{v) = - + CQ + Cxv + 0{v'^), v^Q. (B.9) 

In other words: 

Co - (In^O (-V2|T), - - J^^Z^) - 6 ■ ^ ^ ^ 

Several differential equations can be derived for C,b{'^)- They are all based on the fact 
that c!^^\v) / C,b{v) is 1 and T periodic, thus can be expressed just by matching the 
divergent behavior at the poles with help of Weierstrafi function and its derivatives. 

^^^^^-2 p(.)-p(V2) ^^^^^ (B.ll) 
C^{v) = 2Coa{v) + 2{p{v) - C^)Cb{v) . 

By the same method, one finds a "mirror relation" : 

Cb{v)Cb{-v) = p{b/2) - p{v) (B.12) 

which shows that p{b/2) = Cq — 2Ci. As a consequence, we mention that Cbi^) = 
e~'^°'" Cbi^) satisfies the spin 1 Lame differential equation with spectral parameter 6/2: 

CHiv) - 2p{v)Cb{v) = p{b/2)Cb{v) . (B.13) 
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